Precise charm to strange mass ratio and light quark masses from full lattice QCD 
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By using a single formalism to handle charm, strange and light valence quarks in full lattice QCD 
for the first time, we are able to determine ratios of quark masses to 1%. For m c /m s we obtain 
11.85(16), an order of magnitude more precise than the current PDG average. Combined with 1% 
determinations of the charm quark mass now possible this gives m s (2GeV) = 92.4(1.5) MeV. The 
MILC result for m s /mi = 27.2(3) yields mi(2GeV) = 3.40(7) MeV for the average of u and d quark 
masses. 



Introduction. - The masses of u, d and s quarks are 
some of the least well-known parameters of the Standard 
Model. Even the most inaccurate lepton mass (that of 
the r) is known to better than 0.01% and yet errors on 
light quark masses of 30% are quoted in the Particle Data 
Tables pQ. The reason for the mismatch is the confine- 
ment property of the strong force that obscures the con- 
nection between the properties of the quark constituents 
and the hadron physics that is accessible to experiment. 
To make this connection requires accurate calculations in 
QCD and accurate experimental results for appropriate 
hadronic quantities. A method particularly well-suited 
to this is lattice QCD. Here we will demonstrate its use 
by determining m c /m s to 1% and obtaining as a result 
1.5% errors for light quark masses, which brings them 
almost into line with those of heavy quarks. 

Heavy quark masses, toq, can be determined accu- 
rately because a s (rriQ) is relatively small. 1% errors for 
charm and bottom quark masses have recently become 
possible using 0(a^) calculations in QCD perturbation 
theory for the heavy quark vacuum 'bubble' 2 and there- 
fore for the energy-derivative (or time) moments of cor- 
relation functions for a heavy quark-antiquark pair at 
zero momentum. Since the scale of a s is naturally re- 
lated to the relevant heavy quark mass, the expressions 
can be evaluated accurately. To extract the quark mass 
the perturbative result is compared to a nonperturbative 
determination containing information from experiment. 
For a 1 QQ configuration moments of the experimen- 
tally measured cross-section for (e + e~ — > 7* — > hadrons) 
can be used after isolating the heavy quark contribu- 
tion and using dispersion relations [3] . Alternatively the 
time-moments for heavy quark current-current correla- 
tion functions of various J p can be directly determined 
in lattice QCD calculations that have been tuned so that 
a charmonium or bottomonium mass agrees with experi- 



ment 4, 5 . The time moments must be extrapolated to 
the zero lattice spacing (continuum) limit before the com- 
parison to QCD perturbation theory. These two meth- 
ods give results that agree, with 1% errors for m c (3GeV) 
in the MS scheme. The more traditional 'direct' lat- 
tice QCD method, although somewhat less accurate, also 
gives results in good agreement [5] . We can conclude from 
this that m c is now accurately known. 

The strange quark mass, m s , being much smaller, can- 
not be determined this way and is poorly known at 
present. Instead of a direct determination of m s , how- 
ever, we can use the leverage of an accurate result for the 
ratio m c /m s combined with the accurate m c above [7]. 
But simple ratios of hadron mass differences give unreli- 
able estimates of m c /m s . Two such estimates: 

m(B c ) - m(B u ) = m(S c ) - m(N) = 

m{B s ) - m{B u ) ' m(S) - m(N) ( ' 

differ by almost a factor of 2. The ratio of m s /mi (where 
rn; = (m u + TO,j)/2) is known to about 10% from ra- 
tios of squared masses of K and ir mesons using SU(3) 
chiral perturbation theory [T]. Clearly neither ratio is 
determined well enough this way to provide the accuracy 
we need, because the relationship between hadron mass 
and well-defined running quark mass is more complicated 
than these simple ratios must assume. 

Lattice QCD, on the other hand, can give very accurate 
results for the ratio of two quark masses but only if the 
same formalism is used for both quarks. This has already 
been used to give accurate results for to s /to/, although 
neither m s nor to; is very well determined. Here, for the 
first time, we give an accurate result for m c /m s by using 
the same formalism for charm, strange and light quarks 
and this enables us to cascade the accuracy of the heavy 
quark mass down to the light quarks. 

The Lattice QCD calculation. - Lattice QCD gives 
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direct access to quark masses through the lattice QCD 
Lagrangian. Tuning of the masses is done by calculat- 
ing an appropriate hadron mass and adjusting the quark 
mass until the hadron mass agrees with experiment. Ex- 
perimental measurements of appropriate hadron masses 
are extremely accurate in most cases, with errors at the 
level of tenths or hundredths of a percent. To make maxi- 
mum use of this precision we need to calculate the hadron 
mass in lattice QCD with small statistical and system- 
atic errors. In particular it requires the full effect of sea 
quarks in the hadron to be included. This is now possi- 
ble in lattice QCD [8 . Fixing the four quark masses (mi, 
m s , m c , rrib) from four 'gold-plated' hadrons (n, K, r\ c , 
T) enables other quantities to be calculated with errors 
of a few percent and agreement with experiment is ob- 
tained [HIS]. This is an important test that QCD, with 
only one scale parameter and one mass parameter per 
quark flavor, describes the full range of hadron physics 
consistently. 

The lattice quark mass is a perfectly well-defined run- 
ning quark mass. However, it is scheme-dependent and 
so varies with the discretisation of the Dirac equation 
used in the lattice calculation. For wider applicability 
it is more useful to convert the lattice quark mass to a 
standard continuum scheme such as MS. This renormal- 
ization has been a major source of systematic error in pre- 
vious determinations of light and strange quark masses. 
The best existing result for m s (2GeV), with a 7% error, 
uses the direct method of converting the tuned quark 
mass in the lattice QCD Lagrangian to the MS scheme 
using a 2 s lattice QCD perturbation theory [TO]. The error 
is dominated by the error in the renormalization and it is 
the error that we will remove here, by instead determin- 
ing m c /m s accurately. The Highly Improved Staggered 
Quark action [HJH2] allows us to use the same discretiza- 
tion of QCD for both charm and strange quarks because 
it is a fully relativistic 'light quark' action that can also 
be used for charm quarks. Then the mass renormalisa- 
tion factor cancels in the quark mass ratio. 

We work with eight different ensembles of gluon 
field configurations provided by the MILC collaboration. 
These include the effect of u, d and s sea quarks using the 
improved staggered quark (asqtad) formalism using the 
fourth root 'trick'. This procedure, although 'ugly', ap- 
pears to be a valid discretization of QCD [T5HIT)] . Tests 
include studies of the Dirac operator and comparisons to 
effective field theories. Configurations are available with 
large spatial volumes (> 2.4(fm) 3 ) at multiple values of 
the light sea masses (using m u = m c i = mi) and for a 
wide range of values of the lattice spacing, a. We use 
configurations at five values of a between 0.15 fm and 
0.05 fm with parameters as listed in Table [I] 

On these configurations we have calculated quark prop- 
agators for charm quarks, strange quarks and light quarks 
(again m u = nid = mi) using the HISQ action. The nu- 
merical speed of HISQ means that we have been able to 



Set /3 ri/a auom a ) a l q auam a ) sq L/a T/a N con f x Nt 

1 6.572 2.152(5) 0.0097 0.0484 16 48 631 x 2 

2 6.586 2.138(4) 0.0194 0.0484 16 48 631 x 2 

3 6.76 2.647(3) 0.005 0T05 24 64 678 x 2 

4 6.76 2.618(3) 0.01 0.05 20 64 595 x 2 

5 7.09 3.699(3) 0.0062 0.031 28 96 566 x 4 

6 7.11 3.712(4) 0.0124 0.031 28 96 265 x 4 

7 7.46 5.296(7) 0.0036 (X018 48 144 201 x 2 

8 7.81 7.115(20) 0.0028 0.014 64 192 208 x 2 ~ 



TABLE I: Ensembles (sets) of MILC configurations used, with 
gauge coupling (3, size L 3 x T and sea masses (x tadpole 
parameter Mo) mJJ*'' and rn,Q S s q . The lattice spacing values 
in units of n after 'smoothing' are given in column 3 [14] . 
Column 8 gives the number of configurations and time sources 
per configuration used for calculating correlators. 



Set arrioc 1 + e am Vc amos am Vs 

1 0.81 0.665 2.19381(16) 0.061 0.50490(36) 
0.825 0.656 2.22013(15) 0.066 0.52524(36) 
0.85 0.641 2.26352(15) 0.080 0.57828(34) 

2 0.825 0.656 2.21954(13) 0.066 0.52458(35) 

3 0.65 0.762 1.84578(8) 0.0537 0.43118(18) 

4 0.63 0.774 1.80849(11) 0.0492 0.41436(23) 
0.66 0.756 1.86674(19) 0.0546 0.43654(24) 
0.72 0.72 1.98114(15) 0.05465 0.43675(24) 
0.753 0.70 2.04293(10) 0.06 0.45787(23) 

0.063 0.46937(24) 

5 0.413 0.893 1.28057(7) 0.0337 0.29413(12) 
0.43 0.885 1.31691(7) 0.0358 0.30332(12) 
0.44 0.88 1.33816(7) 0.0366 0.30675(12) 
0.45 0.875 1.35934(7) 0.0382 0.31362(14) 

6 0.427 0.885 1.30731(10) 0.03635 0.30513(20) 

7 0.273 0.951 0.89932(12) 0.0228 0.20621(19) 
0.28 0.949 0.91551(9) 0.024 0.21196(13) 

8 0.195 0.975 0.67119(6) 0.0165 0.15484(14) 
0.018 0.16209(17) 



TABLE II: Results for the masses in lattice units of the 
goldstone pseudoscalars made from valence HISQ charm or 
strange quarks on the different MILC ensembles enumerated 
in Table]!] Columns 2 and 3 give the corresponding bare charm 
quark mass, and Naik coefficient respectively. Column 6 gives 
the bare strange quark mass (e = in that case) . 

use several nearby quark masses for charm and strange 
to allow accurate interpolation to the correct values. Ta- 
ble |H] gives masses for the goldstone pseudoscalar mesons 
made from either a charm quark-antiquark pair or a 
strange one (the r\ c and the T] s ), which are used for tun- 
ing. In the charm case, as well as the quark mass, we 
list the coefficient of the 'Naik' term in the HISQ action 
that corrects for discretisation errors through (amo c ) . 
The quark propagators are generated from random wall 
sources and the goldstone mesons have good signal/noise 
properties so the meson masses can be determined to high 
precision using a standard multi-exponential fit [17] . 

The meson masses can be converted to physical units 
with a determination of the lattice spacing. On an en- 
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1: Grey points show the raw data for every ratio of 
m s on each ensemble (Table In} ; these ratios are fit to 



eq.|4| The dashed line and associated grey error band (and red 
point at a = 0) show our extrapolation of the resulting tuned 
m c /m s to the continuum limit. Blue points with error bars 
are from a simple interpolation, separately for each ensemble, 
to the correct m c /m a , and are shown for illustration. 



semble by ensemble basis this is taken from a parameter 
in the heavy quark potential called r\. Values for r\/a 
determined by the MILC collaboration [H] are given in 
Table [TJ They have errors of 0.3-0.5%. The physical value 
for n must then be obtained by comparing to experimen- 
tally known quantities and we use the value 0.3133(23) 
fm obtained from a set of four such quantities, tested for 
consistency in the continuum limit [T5] [TpJ . 

Using the information about meson masses that we 
have on each ensemble we can interpolate to the cor- 
rect ratio for amo c and amo s using appropriate contin- 
uum values for the masses of the r\ c and r\ s . We cor- 
rect the experimental value of m„ c of 2.9803 GeV to 
w?7 c ,phys = 2.9852(34) GeV. This allows for electromag- 
netic effects (2.4 MeV) [TS] and r/ c annihilation to gluons 
(2.5MeV) [11], both of which are missing from our calcu- 
lation, so increasing the i] c mass. We take a 50% error on 
each of these corrections and also increase the experimen- 
tal error to 3 MeV to allow for the spread of results from 
different rj c production mechanisms [JJ. Since the total 
shift is only around 0.2% of the r] c mass it has a negligible 
effect as can be seen from our error budget below. 

The rj s is not a physical particle in the real world be- 
cause of mixing with other flavor neutral combinations to 
make the r\ and rf . However, in lattice QCD, the particle 
calculated (as here) from only 'connected' quark propag- 
tors does not mix and is a well-defined meson. Its mass 
must be determined by relating its properties to those 
of mesons such as the ir and K that do appear in ex- 
periment. From an analysis of the lattice spacing and 
mj-dependence of the 7t, K, and rj s masses we conclude 
that the value of the rj s mass in the continuum and phys- 
ical mi limits is 0.6858(40) GeV [H]. 

The connection between the MS mass at a scale \i and 



the lattice bare quark mass is given by [TU1 [2D] : 

amo 



m(ji) 



-Z m (na,m a), 



(2) 



Z m = l + a s ( — log(/m) + C + b(amo) 2 

From these two equations it is clear that 

m c {n) am 0c 



am 0s 



(3) 



phys 



where phys denotes extrapolation to the continuum limit 
and physical sea quark mass limit. 

On each ensemble the ratios we have for amo c / am 0s 
then differ from the physical value because of three ef- 
fects: mistuning from the correct physical meson mass; 
finite a effects that need to be extrapolated away and ef- 
fects because the sea light quark masses are not correct. 
We incorporate these into our fitting function: 
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(5) 



are the measures of mistuning, where MC denotes lattice 
values converted to physical units. The last bracket fits 
the finite lattice spacing effects as a power series in even 
powers of a. These can either have a scale set by m c 
(for which we use am Vc /2) or by Aqqd (for which we use 
aro, s ). i,j,k,l all start from zero and are varied in the 
ranges: i, j < 3, k < 6, I < 2 with i + j + k + I < 6. 
Doubling any of the upper limits has negligible effect on 
the final result. The prior on djki is set to 0(1). Sm^t 
is the total difference between the sea-quark masses used 
in the simulation and the correct value for 2mi + m s |18j . 
This has a tiny effect and we simply use a linear term 
(adding higher orders has negligible effect). The prior for 
e? sca is 0.0(1). Figure [l] shows the results of the fit, giving 
m c /m s in the continuum limit as 11.85(16) (x 2 /dof = 
0.42). The error budget is given in Table [TTT| 

m s /mi is known to 1% from lattice QCD as a byprod- 
uct of standard chiral extrapolations of m\ and m 2 K to 
the physical point [5T]. MILC quote 27.2(3) using asq- 
tad quarks [H]. Our HISQ analysis in [T2] gave a re- 
sult in agreement at 27.8(3), using a Bayesian fit to a 
function including terms from chiral perturbation theory 
up to third order in mi and allowing for discretisation 
errors up to and including a 4 and for mixed terms (i.e 
mj-dependent discretisation errors). A full error budget 



is given in Table III the data are given in [TS]. 
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m c /m s 


m s /mi 


overall n uncertainty 


0.4% 


0.1% 


ri/a uncertainties 


0.2 


- 


continuum M Ve 


0.2 


- 


continuum M Vb 


1.1 




Finite volume 




U.o 


a 2 extrapolation, m q interpolns 


0.4 


0.8 


sea-quark mass extrapolation 


0.0 


0.2 


statistical errors 


0.3 


0.4 


Total 


1.3% 


1.0% 



TABLE III: Error budgets for m c /m 3 and m s /mi. 



QCD [HI [53] . These use direct methods of converting 
the lattice mass to the MS mass, and have 10% errors. 
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Conclusions. - Our m c /m s can be used with any value 
for m c to give m s . The best existing result [3] (converted 
from n/=4 to 3) is mi 3) (2.0GeV) = 1.095(11) GeV or 
mi 3) (3.0GeV) = 0.990(10) GeV. Dividing by 11.85(16) 
gives mi 3) (2.0GeV) = 92.4(1.5) MeV and mi 3) (3.0GeV) 
= 83.5(1.4) MeV. 
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FIG. 2: Our results for the 4 lightest quark masses compared 
to the current PDG evaluations (shaded bands) [T]. Each 
mass is quoted in the MS scheme at its conventional scale: 2 
GeV for u, d, s (rif — 3); m c for c (nf = 4). 

Using the MILC values for m s /mi and m u /m,i 
(0.42(4) [14 ) we can then obtain: m| 3) (2.0GeV) 
= 3.40(7) MeV and m| 3) (3.0GeV) = 3.07(6) MeV; 
ml 3) (2.0GeV) = 2.01(14) MeV and m^ 3) (2.0GeV) = 
4.79(16) MeV. The values for all four quark masses are 
plotted in Figure [2] in comparison to the current evalua- 
tions from the Particle Data Tables pQ. 

Thus our high accuracy on m c /m s allows us to lever- 
age 2% accurate values for m s and mi that are com- 
pletely nonperturbative in lattice QCD, for the first time. 
Our m s mass is higher, by around la, than our previous 
value of m s (2GeV) = 87(6) MeV which used 2-loop lat- 
tice QCD perturbation theory [10]. Then the error was 
dominated by unknown a^. terms. Our new result, which 
does not have this limitation, has an error almost five 
times smaller. Our new error is almost an order of mag- 
nitude smaller than other lattice QCD results from full 
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